Vehicular traffic flow at a intersection controlled by signal light with a new 

probability* 
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We introduced a probability of traffic light, Pl, at an intersection when approaching cars in two 
roads are in same conditions. As a application, we proposed a modified Nagel-Schreckenberg cellular 
automata model for describing a confiicting vehicular traffic fiow at the intersection. The results 
show that the plateau region in the fundamental diagrams, caused by the effect of interaction, is 
dependent not only on the probability Pl, but also on the adaptive schemes. 
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Efficient transportation systems are essential for tlie 
every day activities of modern industrialized societies. 
The urgency to reduce CO2 emisions and fuel consump- 
tion, and the excessive, unpredictable travel times dur- 
ing traffic congestion call for more efficient control ap- 
proaches, in particular, the optimization of traffic lights 
in urban road networks. During past years, various traf- 
fic models have been developed to investigate theprob- 
lems of traffic jams [E 0], car accidents 0, 0, U and 
energy dissipation 0] within the framework of single 
lane trafhc models. Recently, a research focus was put 
on urban road networks, which required extending one- 
dimensional traffic models in order to cope with situa- 
tions [1, Q . Evidently the optimization of traffic flow at 
a single intersection is preliminary but important step for 
not only global optimization in city networks but also lo- 
cal optimization. After the first model for simulation two 
crossing roads [lo|, controlling traffic flow at intersections 
have attracted notable attentions 11 1 . 



In principle, the traffic flow at the intersection of two 
roads can be controlled via two distinctive schemes, with 
or without signalized traffic light. It is evident that trafhc 
light is unavoidable as the density of cars increases. And 
it is agreed that a further improvement of the traffic flow 
requires applying flexible strategies than flxed-time con- 
trols. In this letter, we developed a Nagel-Schreckenberg 
cellular automata model for describing a conflicting ve- 
hicular traffic flow at an intersection and present our sim- 
ulation results for three intelligent controlling schemes. 
Our control schemes are inspired by the model proposed 
by Foulaadvand Uj. In that work, approaching car to 



the intersection yield to traffic at the perpendicular di- 
rection by adjusting its velocity to a safe value to avoid 
collision. And the priority is given to the nearest car to 
the intersection according to driving rules. A nature and 
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important question is that which one will get the prior- 
ity when two approaching cars are in same conditions? 
In this letter, we introduce a probability for traffic light, 
Pl, while that happens. In order to capture the basic 
feather of this problem, we have constructed a modified 
NS cellular automata model with three traffic responsive 
schemes applying Pl- 

We first discuss the traffic adaptive controlling scheme 
in which the light signalization is adapted to the traffic. 
Nowadays advanced traffic control systems anticipate the 
traffic approaching intersections. The data obtained via 
traffic detectors installed at the intersection make it pos- 
sible to measure the velocity and distance of an approach- 
ing car and estimate the time the car successively passing 
through the intersection [l^. That supports us to pro- 
mote three possible traffic adaptive controlling schemes. 

In each scheme, the distance to the intersection, the 
velocity and the time passing the crossing point of the 
approaching cars are denoted by di,vi,ti and d2,V2,t2 for 
the first and second road, respectively. Generally, both 
roads are green. By green (red) road we mean the road 
for which the traffic light is green (red). If the two ap- 
proaching cars both can pass the crossing point at the 
next step, one road should change to red to avoid colli- 
sion under the next three algorithms. 

1. The movement priority is given to the nearest 
car((iiniii) to the crossing point, and the car ad- 
justs its velocity as usual with its leading car. In 
contrast, the other road is turned red and the car 
brakes irrespective of its direct gap. The simplest 
way to take into this cautionary braking is to adjust 
the gap with the crossing point itself. If both the 
approaching cars have the same distance, di — d2, 
the priority is given to one road with probability 
Pl, and the other with 1 — Pl; 

2. The priority is given to the car that use less 
time(ti„in) for pass through. While ti ~ t2, the 
priority is given to one road with probability Pl, 
and the other with 1 — Pl. 

3. No matter what the di,Vi,ti is, randomly chose one 
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road and give it the movement priority with prob- 
abihty Pl, and the other with 1 — Pl- 

We do not use the velocity alone to form a new scheme 
not only because of the traffic safety. If we give the pri- 
ority to the fastest car (wmax), the approaching car in 
red road will lost its velocity. For high densities, at the 
intersection, the approaching car in the green road will al- 
ways be faster than that in the red road, and the red road 
will keep in red. With the velocity priority scheme, the 
current of the system rises to its maximum value rapidly, 
and then exhibits linear decrease versus density the same 
manner as in the fundamental diagram of a single road 
(data not shown). 

Referring to Fig. [Tl we consider two perpendicular 
one-dimcnsional chains which represent urban roads ac- 
commodating unidirectional vehicular traffic flow. They 
cross each other at the sites ii = i2=0on the first and 
the second chains respectively. With the closed bound- 
ary conditions , the system is equivalent to cars moving 
in two circles cross at one point with a traffic light. The 
roads, in principle, can each carry two opposite flows of 
vehicles. With no loss of generality, we take the direction 
of traffic flow in the first chain from north to south and 
in the second chain from west to east. Space and time 
are discretizcd in such a way that each chain is divided 
into cells which are the same size as a typical car length. 
Time is assumed to elapse in discrete steps of Is. We 
take the number of cells to be Li — L2 — 500 for both 
roads. Each cell can be either occupied by a car or be 
empty. Moreover, each car can take discrete-valued ve- 
locities 0, I, 2, Wmax in which Umax is the maximum 
velocity of cars. To be more specific, at each step of time, 
the system is characterized by the position and velocity 
configurations of cars and the traffic light state at each 
road. 

The system evolves under the NS dynamics Let 
us briefly explain the NS updating rules which syn- 
chronously evolve the system state from time t to t + 1. 
We denote position, velocity of the ith cat at time step t 
by x{i,t) and v{i,t) , respectively. The same quantities 
for its leading car are correspondingly denoted by -|- 
I, t) and v{i-\-l,t) . The number of empty cells in front of 
the ith car is denoted by d{i, t) — x{i -I- I, i) — x{i, t) — I. 
Concerning the above considerations, the following up- 
dating steps evolve the position and the velocity of each 
car in parallel. 

(1) Acceleration: 

v{i,t+ 1/3) min[v{i,t) + I, Umax]; 

(2) Slowing down: 

v{i,t + 2/3) min[v{i,t + l/3),d{i,t)]; 

(3) Stochastic braking: 

v{i, t-\-l) — > max[v{i, i-|-2/3) — 1, 0] with the probability 

(4) Movement: 

x{i,t+ I) x{i,t) +v{i,t + l). 
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Figure 1: Intersection of two uni-directional roads. They 
intersect each other at the end. A closed boundary condition 
is applied 
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Figure 2: Javg versus density at Pl ~ 0.5 and p — 0. 



The state of the system at time i -I- I is updated from 
that at time t by applying the modified NS dynamical 
rules. After transients, two roads maintain steady-state 
currents, defined as the number of vehicles passing from a 
fixed location per a definite time interval, denoted by Ji , 
J2 and the average Javg — (Ji+ J2)/2. They are functions 
of the global densities pi = Ni/ Li , p2 — N2/L2 and 
P = (pi + P2)/2, where Ni and N2 are the number of 
vehicles in the first and the second road, respectively. We 
kept the global density pi = p2 and Umax = 5 during the 
simulations. Figure 2 exhibits the fundamental diagram 
of the isolate system, Javg versus p, and the difference in 
three schemes at the fixed Pl — 0.5, p = 0. 

It is observed that, in scheme 1 and 3, for small densi- 
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Figure 3: Javg versus p at p 
schemes. 
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Figure 4: Javg versus p at p = for various values of Pl in 
scheme 1. 



ties yO, Javg rises to its maximum value rapidly, and then 
undergoes a short rapid decrease after which a lengthy 
plateau region is formed. The current in the plateau re- 
gion is independent of p. In scheme 2, Javg rises to the 
plateau without any decrease. After the plateau, in all 
three schemes, Javg exhibits linear decrease versus p in 
the same manner as in the fundamental diagram of a sin- 
gle road. It is known that local non-linear interactions 
can lead to system-wide patterns of motion and a lo- 
cal defect can affect the low-dimensional non-equilibrium 
systems on a global scale has been confirmed for cellular 
automata models describing vehicular traffic flow lj| . As 
discussed in the intersection makes the cross point 
as a sitewise dynamical defect whose effect is to form the 
plateau region in our case. In fact, when Pl = 0.5, the 
model yield to a non-signaled intersection and the scheme 
one is equal to that discussed by Foulaadvand. We ob- 
served that the Javg in scheme 2 is higher than others. 
The reason is that, in scheme 2, the car got the prior- 
ity requires less time to pass through. This raises the 
efficiency of the intersection compared to the two other 
schemes hence the second scheme is more optimal than 
them. 

Intersection of two chains makes the intersection point 
appear as a sitewise dynamical defective site, and the 
smaller the Pl is, the stronger the dynamic defect is. 
For each value of Pl, consider the current Ji(J2) in each 
road, the larger the Pl is, the plateau region is wider 
and the current value is more reduceddata not shown. 
In order to find deeper insight, it would be illustrative 
to look at the behavior of average current as a function 
of density. Fig. [3] show that in scheme 3, with fixed 
p — 0, Pl = P2 and < Pl < 0.5. The results show 
that the plateau region of Javg also becomes wider and 
the value becomes smaller as Pl increases. Because of 
the equivalence of the two roads, the condition inverts 



while 0.5 < Pl < 1. In particular, both with Pl = 
and Pl = 1, the fundamental diagram exhibits the same 
manner as that of a single road. The results are same as 
that in scheme 2. But there is little different in scheme 
1. As shown in FiglH with fixed p — 0, pi — p2 and 
< Pl < 0.5, while the length of the plateau region 
grows with the Pl increase, the height of the plateau do 
not show significant dependence on Pl. While for low 
density of roads, Pl should be adapted according the 
situation for higher Javg, for high density, Pl of two roads 
should be set the same, Pl = 0.5, for optimal purpose. 
While for low density, Pl should be adapted according 
the situation for higher Javg 

We then examined Javg versus Pl. Because the two 
roads are equivalent in our model, the symmetry centered 
at Pl = 0.5 is expected, as shown in Fig[5l In general, the 
dependence of Javg on Pl depends on the value of p. For 
large values of p, pci > 0.70 in scheme 1 and Pc2 > 0.61 
in scheme 2 and pcs > 0.66 in scheme 3, Javg increases 
with Pl, then starts its decrease after the maximum at 
Pl = 0.5. For p < pd, the results are different in three 
schemes. In scheme 2, Javg increases with Pl to the 
maximum at Pl2c, then the value reduced untill Pl = 0.5. 
After that point, Javg exhibits an increase up to Pl2c £^nd 
decrease subsequently. We note that the region, Pl G 
[^'l2c, ^L2c]' interested. For each value of p, this region 
is on the same line, though Pl2c and Pl2c different. 
That is the same as shown in scheme 3. In scheme 1, 
there is little different that, for each p except large values, 
Javg is almost independent of Pl and remains constant in 
Pl e [fL2c,fL2c]- In fact, Pl G [Pl2c,Pl2c] denotes the 
appearance of the plateau region as shown in Figl2l The 
results show that, for each p, the height of the plateau 
region is independent on p. In particular, that is also 
independent of Pl in scheme 1. 

As shown in Fig [51 at fixed p = 0.3 and brake possibil- 
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Figure 5: Javg versus Pl at p = 0.3 and p = 



ity p = 0, while Javg changes httle in scheme 1, Javg rises 
to its maximum value up to 0.363 and 0.364 in scheme 
2 and 3 respectively. Then Javg reduces rapidly until 
Pl = 0.5, and the condition inverts subsequently. In 
scheme 1, the most commonly used method in non-signal 
intersections, the value of Javg is less than that in others 
and insensitive to Pl- It is easy to see that the scheme 2 
is the optimal one. While it is difficult for drivers to get 
the velocity and distance to the cross point of another 
approaching car, the traffic light with traffic detectors is 
unavoidable even at Pl = 0.5. 

In summary, we introduced a probability of traffic 
light, Pl, while cars upon reaching the intersection are 
in same conditions. We have investigated the flow char- 



acteristics in a signalized intersection via developing a 
NS model applying Pl. We have considered three types 
of controlling schemes. In particular, we have obtained 
the fundamental diagrams and the dependence of average 
currant on the probability PLroad densities. Our findings 
show the probability of priority at the intersection gives 
rise to formation of plateau regions in the fundamental 
diagrams and an optimal method for controlling the traf- 
fic is unavoidable. The findings may shed light on the way 
to further study the urban cross road. 
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